Abstract 



I introduce a simple continuous probability theory based on the Ginzburg-Landau 
equation that provides for the first time a common analytical basis to relate and describe 
the main features of two seemingly different phenomena of condensed-matter physics, 
namely self-organized criticality and multifractality. Numerical support is given by a 
comparison with reported simulation data. Within the theory the origin of self-organized 
critical phenomena is analysed in terms of a nonlinear singularity spectrum different from 
the typical convex shape due to multifractal measures. 
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The concept of self-organized criticality (SOC) [|T| has attracted great interest re- 
cently, both analytically p, |^ and experimentally §]. The idea behind SOC is that a 
certain class of dynamical many-body systems drive themselves into a statistically sta- 
tionary critical state, with no intrinsic length or time scale, where they exhibit fractal 
behaviour and generate 1// noise. Besides SOC, the generalization of fractal growth 
to self-similar multifractals has also attracted considerable attention from the physics 
community over the past years (see, e.g., P, 0). Theoretical models to describe MF 
have been concerned with mean-field arguments p| and standard renormalization group 
methods 0. These approaches, however, do not allow for the direct analysis of the 
(parabolic-shaped) singularity spectrum typical for a multifractal structure ^ 0. 

Motivated by the recent suggestion that SOC supports the appearence of fractal 
structures it is natural to ask then if there is a common principle underlying the 
seemingly unrelated phenomena of SOC and MF. So far as I know a fixed scale trans- 
formation method |1T|, developed for fractal growth, has been used to investigate ana- 



lytically the nature of 2D clusters in SOC. Henceforth, it is also tempting to search for 
a unifiying scenario that underpins a plausible link between MF and SOC. In fact this 
is the motivation for this work in which I only take a step in that direction. 

In this rapid communication I propose a simple continuous probability theory based 
on the Ginzburg-Landau (GL) equation ]12| that combines together the concepts of 



SOC and MF for the first time. In this goal I explore a new analytical basis which, 
on the one hand, allows to understand the genesis of SOC from the point of view of a 
nonlinear singularity spectrum and, on the other hand, reveals further insight into the 
physics governing this crossover. In short, I make here an attempt to unravel the basic 
phenomena at the origin of power laws and multifractal dimensions. 



A crucial feature of the present formalism is to consider that all random variables 
in a one-dimensional (ID) space, TZ^, are functions of the coordinate variable x which I 
map into an equivalent independent variable ( -say, energy /unit force- characterizing a 
random system. Then, all probabilities may became expressible in terms of the uniform 

probability distribution function \T3 

g{C2) - ^(Ci) = p{Ci < c < C2} ^ 0(c) dc , (1) 

where {} indicates the function interval and is a uniform probability density on the 
line (or 71^) which needs to be specified. 

Assume, within this continuous probability model, that the order parameter satis- 
fies 

0(C) = y{l + /^^(C)} , (2) 

such that 0(C +oo)/0o ~^ and — oo)/0o 1; {i.e., ^H{+oo) = —1 and 

/iif(+oo) = 1), and is a coefficient as discussed below (c./., Eq.([T9D). 

I postulate H{Q in Eq.(^ to be given by the real solutions of the following static, 
dimensionless GL-like equation 0: d'^H{C)ld{Cf + pH{C) - qH^iC) = 0, C*^^, 
where (* = C/Coj D is the ^-domain, [p,q] > and Co is a positive coefficient of 
dim[length]. It is not unreasonable to consider the (ID) domain D to be infinite. Then, 
the GL posseses the well-known stable, kink solution 

H{C) = ±.l^tanh{C\ll) . (3) 

y 5 V z 

Using this result together with Eqs.(|I]) and (0), I will establish a simple relation for the 
probability distribution P. 
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The integral of Eq.(|l|) over the hmits: C2 = ^iCo > C + -^2^ = Ci is 

^(AiCo) - ^(C + A2C0) = ^ /''^° {1 ± /X tanh£)} dC ^ -r(C) , (4) 

where I introduced the new function t{Q and set p = g = 2 to reduce the number of 
free parameters. These integration hmits lead to the condition 

A2-Ai + C*<0 . (5) 

Note that the sign in Eq.@ implies that the single functions Q are here assumed to 
satisfy the condition ^(C + A2C0) > ^(AiCo) for C 7^ 0, which throughout the theory are 
undefined, whereas the free parameters Aj (z=l,2) will restrict the range of via Eq.(^). 
Suppose A2 7^ Ai, then it is straightforward to solve the above integral and get 

r(C) ^ (1 + -^— ){r(O)T/i0*lncos/iA2} 
A2 — Ai 

^f''4>o{^ — Incos/iAi + lncos/i(A2 + C*)} 5 (6) 

A2 — Ai 

in which 0q = Co0o/2 and 

r{0)^g{X,Co)-giXiCo) = iX2-XiWo{lTf^rx} , (7) 
such that 

^ Incos/iAi — lncos/iA2 , ^ 

Ta^ ■ (8) 

A2 — Ai 

When (* = 0, MF will be described via Eq.(^ assuming ^(A2Co) < ^(AiCo), whereas for 

soc: e;(A2Co) > QiXiCo). 

To analyse the possible multifractal features of P let me define the function D^* as 
r{C)^{X,-X^ + C}Dc* . (9) 
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From Eqs.@ and (Q) it follows that 

-D^* = T — ^— —{t (0) ^ ^(pl In cosh X2} 
A2 — Ai 



± , ^^ -^\ncosh\^ + In cosh{X, + C)} , (10) 

A2 — Ai + C, A2 — Ai 



such that A2 — Ai + C* 7^ 0. From this relation I obtain 
r(0) 



A2 — M 



and 



De^_oo = /^c*-oT/i0S{i-rA} . (12) 

However if A2 — Ai + C* = then 

-Da, :( -1 1 = ± Ai^ol^A + ton/iAi} . (13) 

Complementary to Eq.(^ I also define the following dimensionless function 

a{C) = ^^(C*) = ^c*-o ± /i0S{rA + tanh{X2 + C)} ■ (14) 
Therefore, using Eq.(|l^) in conjunction with Eq.(^), it can be easily shown that 

"max = a(C* -00) = D)^*-.-oo , 

"mm = a(C* +00) = D)^*^+oo • (15) 

In order to relate these equations to the principle of MF |T1|, I consider the following 
condition for a and D. If (* — ^ +00, then -D^^+oo and «(C*) are allow to take one of 
the two values, namely or 2L)^.^o, depending on the phenomena to be considered. The 
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particular case for which a{(* —>■ +00) ^ and -Df*^+oo ~ will be in correspondence 
with the features of MF. To achieve this I approximate 

where the integer factor e = ±1 is included to distinguish these two limiting values. 



According to the original definitions used in MF [|T5], [T6| I also define the dimension- 
less functions 

/(a)^C«(n-r(n , (17) 

and 

Cc^-Q^riC) = Tf^<P*oSech\\2 + n • (18) 

When f{a) and D(^* are smooth functions of a and (*, then /(a) can be related to t{(*) 
by a Legendre transformation ||T0[ . As is well known this reflects a deep connection with 
the thermodynamic formalism of equilibrium statistical mechanics where r and C(^* are 
the analogous 'free energy' and 'specific heat', respectively. Within this characterization 



a becomes an analogous 'internal energy' and / an analogous 'entropy' [|15 . 

Having established a physical meaning for these functions and assuming 0q > as 
discussed below, I realize then that fi can take the values 

/i^Tl , (19) 

since it is physically satisfactory to set C^. > 0. This, in turn, implies that the reduced 
probability density 0(C) /0o of Eq.(0) can take the desired values and 1 when (* — > +00 
and C* —00, respectively. Within these limits, in the case of MF, the quantities a{C*) 
of Eq.dg), D^, of Eq.([Tg) and of Eq.(|g) converge to zero whereas r(C*) of Eq.(§ 
increases or decreases monotonically dependending on (*. 



Besides this, the choice of Eq.(|19D imphes that I may also obtain Eq.(^) from Eq.(|TT]) 
in conjuction with Eq.(^ provided that e +1. However, in the case that e ^ —1 
is choosen, Eq.(|16[) can only be recovered using again the aformentioned equations but 
changing the sign of the previous value for t(0). This will became more clear latter 
in the calculations. The present continuous probability theory is thus dependent on A, 
(2=1,2) and (the sign and magnitude of) r(0) in Eq.(|^), where the reduced variable C, 
satisfies the condition given by Eq.(|^) and 0o is positive satifying Eq.([T6|). 

The case Ai > A2 such that A2 > (MP): In Fig. 1(a) I display the dependence of the 
analogous 'free energy' r on the reduced coordinate variable C,* for different values of Ai 
and A2, such that r(0) = —1 and e = 1. Noting that the difference between A2 and Ai 
is —1 in all three curves illustrated, so as to have -D(^*_>o = 1? then from Eq.(^) C ^ 1- 
Hitherto, the present GL-based approch to MP allows C* to take on negative as well as 
positive values < 1. However, in the following plots concerning MP, the range of C,* is 
extended to 3 for illustrative purposes. 

Prom Pig. 1(a) it can be seen that, on increasing the value of Ai, there is a more 
rapid convergence of t((^*) for positive C,* than within the region C* < 0; displaying 
thus typical features of MP |T5|. Such a nontrivial behaviour of r illustrates the data 
collapse or breakdown of MP at C* > where t{C*) > 0. This is in accordance with 
the MP structure of the function D^. shown in Pig.2. In fact, this corresponds to a 
spectrum of fractal dimensions. In this figure I point out that D^* -defined via Eq. (p!oD - 
is not constant for positive and negative values of (* so that it may correspond to 
a multifractal dimension. Hence, r becomes also a nonlinear function of (*. In fact 
this approach yields, e.g, for the case Ai = 1,A2 = (full line in Pig. 1(a)) the values 
D^*^i = 0.421 and Z)^.^_oo = 3.53 using Eqs.([T3|) and (p!2D, respectively. 

6 



On the other hand, in Fig.3(a) I display the analogous 'internal energy' a as a 
function of (* for the same set of values of Ai and A2 as in Fig. 1(a). When Ai = 1 and 
A2 = this function exhibits sharp variations around = with a maximum value 
that shows a stronger dependence on Ai for positive (* than for negative (*. This rules 
out the possibility that the actual position for a critical value of ( (or analogous critical 
inverse 'temperature' ||T3|), at which the multifractal formalism actually breaks down, is 



obtained at (* = 0. For other values of Ai and A2, C changes towards negative values. 
One should also bear in mind that for (* < (*, the analogous 'free energy' r(C*) is 
dominated by a which, in turn, varies with the magnitude of A2. 

Characteristic features of a phase transition at Q can be figure out by examining 
the shape of the analogous 'specific heat' C^* of Eq.(p!8D which is illustrated in Fig.4. 
It can be easily visualized that there is a sharp peak around the value Q = for the 
case corresponding to the full line in Fig. 1(a). The heights and positions of these curves 
are strongly dependent on Ai > A2. Nicely, this finding is also similar to reported MF 
results [0. Finally, I investigate the behaviour of the analogous 'entropy' / against a 



for several values of Ai and A2. As can be seen from Fig. 5(a) I find that, on incresing the 
magnitude of Ai, the left-hand side of these plots converge more rapidly than the right- 
hand sides which converge poorly. This is in complete agreement with the MF signal 
observed in the context of self-similar random resistor networks (open and full circles 
in Fig.5(a) from |0) or, e.g., in DLA when Ai is taken to be related to the system 
size of the simulation box 0, 0. In fact, this singularity /(a) shows the characteristic 



normalized convex shape found in MF. Given its interpretation of a multifractal character 
for fractal subsets, each with a different fractal dimension having singularity strength a, 
it is expected that /(a) > 0. 



Let me remark again that all of these predictions resemble qualitatively many of the 
intriguing results observed in MF, such as the breakdown of multifractal behaviour and 
the existence of a phase transition at C*, which may not necessarily be at C* = 0. In the 
examples considered above the maximum and minimum values of the analogous 'internal 
energy' a can be estimated from Eqs.(P3D in conjuction with Eq.(|T^, respectively. This 
allows for the existence of a critical point C* above which the infinite hierarchy (or 
broad distribution) of phases can be found, but below which a single phase appears 
characterized by the maximum 'energy' amax- 

The case Ai > A2 such that A2 < (SOC): In Fig. 1(b) the dependence of the normal- 
ized probability distribution function r of Eq.(|^) on the energy /unit force variable C,* 
is plotted for different values of Ai and A2 such that, as a main difference with respect 
to MF, r(0) = 1 and e = —1. This new choice of r(0) allows to normalize t{(*) and the 
choice of e enables to mimic the main features of SOC, namely a power-law behaviour, 
provided the reduced variable (* is associated with the /o^f-function of some measured 
event. As an illustrative example, in Fig. 1(c) I show a reasonable description of the 
probability vr that the water flow intensity of a randomly chosen site is larger than s, 
as obtained from the SOC signal calculated within a model of erosion |T^. This de- 
mostrates that the theory do apply to self-organizing systems. For a class of continuous. 



cellular automaton models of earthquakes [|19| , |20[| , the function r is simply reinterpreted 
as being the number of events with reduced released energy ~ e^* . 

The several theoretical curves in Fig. 1(b) refer to different values of ( ) Ai = 3 

and A2 = -6; ( ) Ai = 3, A2 = -7; ( ) Ai = 3, A2 = -8. Using Eq.(|), 

then the GL-based theory is valid for C ^ 9, 10, 11, respectively, where I deal with the 
physically interesting interval < t{(*) < 1. A given slope of the linear behaviour of 
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the curves in this figure is determined by fixing Ai > 0. In particular, this parameter 
may be associated with the elastic parameter of the spring-block model for earthquakes 
which links the rate of occurrence of earthquakes of magnitude M greater than m to 
the energy (seismic moment) E released during the earthquake via the famous empirical 
Gutenber-Richter law ||T^. And the cutoff in the axes C,* may be related to the system 
size of cellular automaton modeling for arrays of threshold elements. 

To see more clearly power-law features in the behaviour of r over a wide range of 
positive values of C,* (as those in Figs. 1(b) and (c)), it is necessary to investigate the 
behavior of the derivative of t(C*), defined through a of Eq. ([I^) , which is plotted in 
Fig. 3(b). A glance into the behaviour of «(C*) indicates that for the smallest displayed 
values of C* > the a-function converges to a constant negative value, revealing in 
this way the constant nature of the negative slopes in the (t-^*) curves of Fig. 1(b). 
On increasing C* each curve smoothly approach a smaller value. Clearly, due to the 
probabilistic nature of r for SOC {i) such convergences of a towards smaller negative 
values need not to be considered and (u) the relations between amax^min and -D(*^ipoo 
(c./., Eq.(p!5|)), become meaningless. 

After establishing this resemblance of a power-law description, it is tempting to con- 
tinue applying anew the above formalism of MF to analyse SOC. As already mentioned, 
this may help to understand the origin of SOC from the novel viewpoint of a nonlin- 
ear singularity spectrum different from what is common to multifractal objects (c./., 
Fig.5(a)). 

In view of the features in q;(C*) of Fig. 3(b), the second derivative of r (not shown), 
namely C^* of Eq. (p^Sf) , present a sharp peak around an inflection point of the function 
a{C,*) of Fig.3(b), say Q^j. As a difference to the above results of MF (c./.. Fig. 4), 
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these peak heights reduce their magnitude on decresing A2 < Ai and shift their position 
towards positive values of (*. In this case no phase transition as in the case of MF is 
expected because r > restricts the range of vahd (* < Q^j. Moreover, as mentioned 
above, D^* remains undefined for SOC due to the re-interpretation of t((^*). Similarly 
to what is shown in Fig. 2 this function for SOC may simply become not constant on 
increasing C^* as in MF. 

On the other hand, the continuous singularity spectrum f{a) plays an alternative 
role when dealing with SOC which can be assessed from Fig.5(b). In this plot it can be 
seen that this function exhibits a rather nonlinear behaviour different from the parabolic 
behaviour typical of multifractal entities (c./., Fig.5(a)). In MF / takes its maximum at 
the value a(C* = 0) whereas in SOC this spectrum becames a monotonically increasing 
(negative) function of (negative) a. On decresing the magnitude of A2 these curves 
converge to —1 and separate out as a function of decreasing a whithin the range of 
validity of a in Fig.3(b). Since r(^*) is positive, then —2 < /(a) < —1. I suggest this 
new aspect of the behaviour of /(a) to be a fundamental property for the additional 
characterization of SOC. It is therefore most likely that the linear behaviour displayed by 
r(C*) (c./.. Fig. 1(b)), that is quantified via a(C*) (c/.. Fig. 3(b)), finds its root through 
the behaviour of Eq.([r^) for /(a) (c./.. Fig. 5(b)). 

To some extent, I have been able to shed light on a basic mechanism leading to 
both concepts of MF and SOC from the continuous density probability as given 
in Eq.(0). This function has been assumed to be related to H{Q via Eq.(H), which I 
postulated to be given by the real kink solutions of a dimensionless GL-like equation. 
Such solutions are known to minimize the GL free energy functional when | p is related 



to a particle concentration [|T^. Of course to use a continuous probability theory may 
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be seen as being heuristic, specially so if reported simulation results have been done 
using discretized cell configurations. But, as I have discussed, a great deal of relevant 
information can be extracted from a continuos approach which, essentially, does relay 
on the sign of A2, t(0) and e only. 

The thermodynamic treatment of MF provides a convenient way to quantify the 
complexity of multifractals by characterizing scale-invariant singularities and, in most 
cases, allowing for an effective comparison between theory and experiments |0, |14|, [16 



Nevertheless, I have proposed here that this formalism may be extended as a tool to 
rationalize SOC on the basis of a singularity spectrum that turns out to be nonlinear as 
well (when comparing Figs. 5 (a) and (b)). While the present (static) GL-based theory 
is extremely simple, it is important to emphasize that it gives information about the 
complex origin of self-organized critical phenomena whose physics is assumed -to a good 
approximation- to be analog to that required to describe MF. This theory also reflects 
the minimal ingredients than can give rise to an intrinsically critical state. In passing, I 
add that multifractal structures can also emerge from GL equations with random initial 
conditions for the temporal evolution [pT| . 

It is a pleasure to acknowledge helpful remarks by Drs. V.L. Nguyen, G.C. Barker, 
W. Wang and H. Rosu. 
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ure captions 

Figure 1: (a) Dependence of the analogous 'free energy' r on the coordinate 
variable (* displaying typical features of MF phenomena Different values of 

( ) Ai = 1 and A2 = 0; ( ) Ai = 1.5, A2=0.5; ( ) Ai = 2, A2 = 1, 

illustrates the data collapse at > where r > 0. In these examples r(0) = — 1 
and e = 1. 

(b) Dependence of the normalized probability distribution function r on the vari- 
able C* displaying typical features of SOC. Different curves refer to different values 

of ( ) Ai = 3 and A2 = —6; ( ) Ai = 3, A2 = —7; and ( ) Ai = 3, 

A2 = —8. In these examples r(0) = 1 and e = —1. 

(c) Present description of the probability vr that the water flow intensity of a 
randomly chosen site is larger than the event s. Open circles are the SOC signal 
in a model of erosion [p!8| . 

Figure 2: Plot of the multifractal dimension D^* against (* showing a non- 
constant behaviour for the same set of values for Ai and A2 as in Fig. 1(a). 

Figure 3: (a) Dependence of the analogous 'internal energy' a on (* displaying 
characteristic features of a phase transition for the same set of values for Ai and 
A2 as in Fig. 1(a). 

(b) The function a of Eq. (p!^ as a function of (* for the same set of values for Ai 
and A2 as in Fig. 1(b). These curves reflect the power-law features in the behaviour 
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of t{C) for SOC. 



• Figure 4: Analogous 'specific heat' C^* as a function of (* displaying characteristic 
features of a phase transition for the same set of values for Ai and A2 as in Fig. 1(a). 

• Figure 5: (a) Dependence of the analogous 'entropy' / on the analogous 'internal 
energy' a for the same set of values for Ai and A2 as in Fig. 1(a). Open and full 
circles are the MF signal calculated in the context of self-similar random resistor 
networks for different cell sizes (see ||17|| ). 

(b) The nonhnear /(a) singularity spectrum obtained from Eq. ([T7|) proposed to 
characterize SOC using the same set of values for Ai and A2 as in Fig. 1(b). 
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